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ABSTRACT 

We present the results of the SuperLupus Survey for transiting hot Jupiter 
planets, which monitored a single Galactic disk field spanning 0.66 deg 2 for 108 
nights over three years. Ten candidates were detected: one is a transiting planet, 
two remain candidates, and seven have been subsequently identified as false pos- 
itives. We construct a new image quality metric, Sj, based on the behaviour of 
26,859 light curves, which allows us to discard poor images in an objective and 
quantitative manner. Furthermore, in some cases we are able to identify statis- 
tical false positives by analysing temporal correlations between Sj and transit 
signatures. We use Monte Carlo simulations to measure our detection efficiency 
by injecting artificial transits onto real light curves and applying identical se- 
lection criteria as used in our survey. We find at 90% confidence level that 
0. 10^0^3% of dwarf stars host a hot Jupiter with a period of 1-10 days. Our re- 
sults are consistent with other transit surveys, but appear consistently lower than 
the hot Jupiter frequencies reported from radial velocity surveys, a difference we 
attribute, at least in part, to the difference in stellar populations probed. In 
light of our determination of the frequency of hot Jupiters in Galactic field stars, 
previous null results for transiting planets in open cluster and globular cluster 
surveys no longer appear anomalously low. 

Subject headings: Planets and satellites: detection - Techniques: photometric 
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1. Introduction 

The discovery of short-period, giant extrasolar planets (Mayor & Queloz 1995; Marcy 
& Butler 1996) provided the exciting potential for large numbers of planets to be discov- 
ered by the transit method, as these "hot Jupiters" have a ~10% geometric probability of 
transiting, and do so every few days. However early predictions greatly overestimated the 
actual discovery rate (Home 2003). This discrepancy resulted from simplistic assumptions 
and a misunderstanding of the effects that systematic noise would play in lowering detec- 
tion efficiency (Pont et al. 2006). It also resulted from an over-estimation of the frequency 
of hot Jupiters. More accurate, and lower, predictions were provided in Beatty & Gaudi 
(2008), where it was noted that objective and quantifiable detection criteria were required 
for more robust inferences of planet frequencies. By adopting such objective and quantifi- 
able detection criteria, we propose that the frequency of hot Jupiters in the field has been 
overestimated by a factor of three (Bayliss & Sackett 2011). 

The SuperLupus Survey was established to detect hot Jupiters in a field positioned just 
above the Galactic plane (b = 11°), and also to determine the fraction of stars that host hot 
Jupiters in this typical Galactic field. The survey design and data analysis was constructed 
so as to fulfill both of these objectives. 

In Section 2 of this paper, we describe the SuperLupus observations and data reduction. 
The photometry is described in Section 3. The criteria for candidate selection are detailed 
in Section 4, and the 10 identified candidates are described and analyzed in Section 5. 

In Section 6, we set out the details of the Monte Carlo simulations performed to calculate 
the detection efficiency and the effective number of stars probed for planets in the SuperLupus 
Survey. In Section 7, we apply this efficiency to the actual results of the survey to determine 
the fraction of stars in the field that host a hot Jupiter. Finally, in Section 8 we summarize 
and discuss the implications of our results. 

2. Observations and Data Reduction 

The SuperLupus Survey is an extension of the Lupus Survey (Bayliss et al. 2009), in 
which the duration of the original survey was approximately doubled to provide sensitivity 
to transiting planets with periods as long as 10 days. Many details of the data reduction 
procedure are set out elsewhere (Bayliss et al. 2009); here we summarize the procedure and 
highlight aspects that differed between the Lupus and SuperLupus projects. 

Both surveys monitored a single field using the Wide-Field Imager on the ANU 40- 
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Inch Telescope at Siding Spring Observatory. The field, in the constellation of Lupus, is 
11° above the Galactic plane centered at R.A.=15 h 30 m 36.3 s , Decl.=-42°53'53.0" (J2000). 
In total, 5158 images of 300 s exposure were taken in a custom (V+R) filter. Of these, 
2801 images were from the original Lupus Survey (2005 and 2006) and an additional 2357 
new images were taken in 2008. The average full-width half-maximum of the point spread 
function (PSF) over all 5158 images is 2.02". 

Initial data reduction (bias, dark, and flatfield correction) was carried out on the mosaic 
frames using standard IRAF 1 tasks in the package MSCRED. Sky flats were obtained at 
twilight whenever conditions permitted, and master flatfields were produced by median com- 
bining the ~30 sky flats most proximate in time to each night. Bad pixels and columns were 
masked and the mosaic frames were split into individual CCD frames. A world coordinate 
system solution was calculated for each image in order to identify stars for photometry. 

3. Photometry and Detrending 

Photometry for the original Lupus Survey was performed using Difference Imaging Anal- 
ysis (DIA: Alard & Lupton 1998; Wozniak 2000). For the SuperLupus Survey, we instead 
used Source Extractor (Bertin & Arnouts 1996) aperture photometry. This approach was 
motivated by the fact that many of the survey images have PSFs that are asymmetric due 
to tracking or focusing issues. Such PSFs can cause potential problems for DIA photometry, 
but are not of concern for aperture photometry as long as the photometric aperture is large 
enough. Initial tests revealed that aperture photometry could achieve the same (and even 
better) photometric precision for the majority of the brighter stars (V < 20) in the field. In 
fact we note that the final SuperLupus aperture photometry resulted in approximately 10% 
more stars with a root-mean-squared (RMS) variability less than 0.025 mag compared with 
the original Lupus Survey using DIA. 

We found that a photometric aperture radius of 10 pixels (3.75") resulted in the highest 
precision lightcurves for the majority of these stars. 

A reference catalog of 50,907 target stars (V < 20) was produced from a single, high 
quality, reference image. Aperture photometry was performed on these stars for each survey 
image, with apertures re-centered on the point of peak flux for each star in each image. 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Background subtraction was performed using a background map produced for each image 
by median filtering and bi-cubic-spline interpolation over a 8x8 pixel grid. 

Systematic trends in the resulting lightcurves were removed using the Sys-Rem algorithm 
(Tamuz et al. 2005) with 12 iterations. To further refine this dataset, we removed stars with 
lightcurves that displayed RMS variability greater than 0.05 mag. We also discarded poor 
quality images as described below. In total this left us with 26,859 stars over 3585 images. 

3.1. Image Quality Metric 

To identify poor quality images, we define an image quality metric (Sj) as: 

i n 

where for the i th star on the j th image, is the average subtracted residual magnitude 
and Oij is the photometric uncertainty. This metric is computed in the Sys-Rem algorithm, 
where linear systematic effects are removed in order to minimise Sj (Tamuz et al. 2005). In 
essence, Sj is a measure of how all the stars in image j vary from the mean of all images 
in the survey. We retained only those images with Sj < 0.02 for searching for planetary 
transits. The advantage of using this metric is that it uses information from all stars over 
the chip in a single quality factor that is directly related to the precision of the resulting 
stellar photometry. 

This image quality metric also proved useful to determine if transit candidates were 
statistical false positives by looking for correlations between "transit" events and Sj (see 
Section 5). 

4. Transit Search and Detection Criteria 

In order to accurately determine the efficiency of the SuperLupus Survey, candidates 
must be identified in an automated and systematic manner that can be applied identically to 
synthetic lightcurves in the Monte Carlo simulation (see Section 6). This was implemented 
by way of a set of six selection criteria that were automatically applied to each processed 
light curve. 

Initially all 26,859 SuperLupus lightcurves were searched for transit events using the 
Box-fitting Least-Squares (BLS) algorithm (Kovacs et al. 2002). Trial periods in the range 
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1.01 < P < 10 d were tested with 55,000 equally-spaced frequency steps and 200 phase bins 
per trial frequency. Candidates were then identified based on their BLS Signal Detection 
Efficiency (SDE), as defined in Kovacs et al. (2002). After testing synthetic lightcurves 
with the same temporal and noise characteristics as our data, a candidate threshold value 
of SDE = 7.0 was adopted. 

We defined the effective signal-to-noise of the detected transit event, a, for all candi- 
dates, as: 

a = -y/n^Ts , (2) 

(7 

where 5 is the BLS-determined transit depth, a is the RMS variability of the entire (Sys-Rem 
corrected) lightcurve, and n obs is the actual number of data points within the BLS-determined 
transit event, i.e. the number of "in-transit" data points. A threshold of a = 10 was adopted, 
identical to the threshold adopted in two other transit surveys of faint stars (Burke et al. 
2006; Hartman et al. 2009). At this stage we also rejected candidates for which more than 
90% of in-transit data originated from an event on a single night. 

As with other ground-based surveys, many of our candidates displayed integer-day pe- 
riods resulting from nightly systematic effects. Therefore candidates with 1.01 < P < 1.02 
and 1.98 < P < 2.02 were rejected. We also required that candidates had a magnitude of 
V < 18.8. 

Finally candidates were cross-matched with the catalog of known variable stars in the 
field (Weldrake & Bayliss 2008) and eclipsing binary systems were removed. 

The selection criteria are summarized in Table 1. 



5. Initial Candidates 

A total of 10 lightcurves fulfilled the six selection criteria set out in Section 4. Their 
characteristics are described in Table 2. 

The candidate SL-07 was immediately recognised as Lupus-TR-3, a star with a transiting 
planet that had been discovered during the original Lupus Survey (Weldrake et al. 2008; 
Bayliss et al. 2009). To investigate the likely nature of the other candidates, we used the r) p 
diagnostic method proposed by Tingley & Sackett (2005). The results are set out in Table 3. 
Candidates SL-02 and SL-04 show very high diagnostic numbers, indicating that their transit 
parameters make them unlikely to be real transiting planets. In the case of candidate SL-04, 
a very clear transit was observed on 6 June 2005, with a transit time of approximately 6 hr, 
in contrast to the BLS determined duration of 8.08 hr. Use of the shorter duration reduces 
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the r) p diagnostic value to 1.9, which, while still high, could be explained by a large radius 
planet, since the rj p diagnostic is normalized to the assumption R P = lRj. 

To help identify statistical false positives, we checked if the detected transit event was 
correlated with the image quality, using the calculated value of the metric Sj, as defined in 
Eq. (1). 

Three candidates, SL-08, SL-09 and SL-10, were found to have transit events that were 
strongly correlated with Sj, indicating that the detected transit signature is almost certainly 
a systematic effect that the Sys-Rem algorithm had not fully removed. An example of such 
a correlation is shown in Figure 1. 

We inspected the images of the candidates for signs that the transit signature might 
be caused by proximity to the edge of a CCD, a bad column, or a bright star. We also 
inspected the image for evidence of significant flux contributions from neighboring stars 
within the photometric aperture. 

SL-01 was found to be blended with a very bright, saturated neighboring star at a 
distance of 27". We therefore ruled this out as being a genuine candidate. We found that 
the candidate SL-03 had three neighbors close to the photometric aperture, and that SL-05 
had a neighbor well within the photometric aperture. Photometry for these candidates was 
re-derived using a smaller photometric aperture to avoid the contaminating neighbors. This 
analysis revealed that the transit event seen in SL-05 was in fact due to a very deep eclipse 
(40%) occurring in a faint neighbour situated just 2.6" from SL-05. It also revealed that 
while the SL-03 target star was responsible for the observed transit event, close neighbors 
and the large aperture of the original photometry had diluted the size of the dip from the 
true depth of 0.100 mag to a shallower 0.042 mag. The true depth is characteristic of an 
eclipsing binary and does not pass our selection criterion. 

In summary, we conclude that of the 10 candidates identified in the SuperLupus Survey, 
only SL-07 can be confirmed as a genuine transiting planet. We rule out seven candidates 
(SL-01, SL-02, SL-03, SL-05, SL-08, SL-09, SL-10) as being false positives, while two candi- 
dates (SL-04 and SL-06) remain. These two remaining candidates are V=17.8 and V=18.6 
respectively, making follow-up extremely difficult. Consequently, we have not observed these 
candidates further, although formally we cannot exclude them. In Section 7 we calculate 
hot Jupiter frequencies assuming only SL-07 is a transiting planet, but also provide figures 
for the case where one of these remaining candidates is a transiting planet. 
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Fig. 1. — Phase- wrapped lightcurve for candidate SL-08 (top panel) and the corresponding 
image quality, Sj, for images with Sj > 0.005 (bottom panel) wrapped at the same period 
and phase. The correlation between the transit event and image quality indicates that SL-08 
is a statistical false positive. 
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6. Monte Carlo Simulations to Determine Detection Efficiency 

To determine the efficiency of the SuperLupus Survey to detecting transiting planets, we 
carried out a Monte Carlo simulation of the SuperLupus Survey. Synthetic transit signatures 
were inserted into actual SuperLupus lightcurves and then recovery attempted using the same 
methodology and selection criteria used for our survey. 



Generating a realistic set of synthetic transit lightcurves against which to test a detection 
algorithm requires knowledge of the actual distribution of transit depths, durations, and 
periods of transiting planets. Ideally, this would be achieved by determining the radius of 
each monitored star in the field. This is possible when monitoring an equidistant population 
such as an open cluster or a globular cluster (e.g. Hartman et al. 2009; Weldrake et al. 2005). 
Since the distances to field stars cannot be determined from photometry alone, the exact 
radius of each star in the SuperLupus field is unknown. We therefore used the Beasongon 
model of the Galaxy (Robin et al. 2003) to provide a statistical distribution of stellar radii 
and masses for the stars monitored in our field. The population was synthesized over a 
0.66 deg 2 field centered at b = 11°, / = 331°. Interstellar extinction was set at 0.6 magkpc -1 , 
based on the Schlegel dust maps of the region (Schlegel et al. 1998). Stars were selected in 
the magnitude range 14 < V < 18.8 to match the SuperLupus observational and selection 
constraints. The Beasongon model returned a population of stars that closely matched both 
the total star count, and the distribution of apparent magnitudes (see Figure 2), giving us 
confidence this model accurately simulates our field. 

We note that the Beasongon model indicated that 24% of the survey stars have log g < 
4.0. In our Monte Carlo these stars are classified as giants and given zero probability for 
transit detection. We assign a radius to each star in the simulation based on the stellar mass 
from the Beasangon model and the simple relationship given by Cox (2000): 



6.1. Modeling the Stellar Population 



log— = 0.917 log 



M 



- 0.020 . 



(3) 



M 



© 



6.2. Lightcurve Generation 



Lightcurves for the Monte Carlo were created by first randomly selecting a star from 
the Beasangon model population. The temporal sampling of the star is set to the actual 
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Fig. 2. — Comparison of the apparent V magnitude distribution between the Beasoncon 
model (dashed line) and our observed field (solid line). 

timestamps of the SuperLupus Survey. We then need to attach an appropriate photometric 
uncertainty to each point on the synthetic curve. 

This was done by using the V magnitude of the chosen star to randomly select an RMS 
uncertainty from the actual distributions in RMS observed in the SuperLupus Survey for 
stars of that V magnitude, as described below. 

The RMS scatter for each real lightcurve in the SuperLupus Survey, after the application 
of Sys-Rem, shows a wide range in RMS scatter at a given V magnitude. Therefore we 
calculated the actual distribution of RMS scatter for each V magnitude in the survey, as 
displayed in Figure 3. As an example, the distribution for V = 17.72 is shown in Figure 4. 
It can be thought of as a vertical slice through the density plot in Figure 3. 

It is evident that the RMS scatter distributions are non-Gaussian, with long tails to- 
wards high RMS scatter. This further highlights the importance of basing the RMS scatter 
of synthetic lightcurves on the actual distributions of photometric precision rather than ap- 
proximating them using a single-valued function of magnitude as has been done in some 
previous studies (e.g. Weldrake et al. 2005; Gould et al. 2006). 

A value for the photometric precision is selected randomly, using the appropriate RMS 
scatter distribution as a probability weighting, for each Beasongon star. Once the value of 
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Fig. 3. — The number density of stars with a given RMS photometric precision plotted 
against V mag. The number density is given on a logarithmic scale. The bin width is 
0.04 mag for V magnitude, and 0.3 mmag for RMS photometric precision. 



- 11 - 




Fig. 4. — RMS photometric precision distribution profile for stars in the SuperLupus Survey 
between V = 17.70 and V = 17.74. 
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the photometric precision is determined, a real SuperLupus lightcurve is selected with the 
same photometric precision, and used as the base lightcurve into which we inject a transit. 

6.3. Injecting Transits 

To simulate the lightcurve from a transiting hot Jupiter, we first determined the param- 
eters for each simulated system. The stellar mass and radius were taken from the Beasoncon 
model and Eq. 3. The planetary mass was fixed at 1 Mj, but we note that this does not 
significantly alter the transit parameters when Mp << M*. Four different planetary radii 
were simulated: Rp = 0.8, 1.0, 1.2, and 1.4 Rj. 

The inclination of the orbit, i, was determined by choosing a random cosz between zero 
and cosi min , the limiting inclination that results in a transit (Sackett 1999): 

R* + Rp 

COS l m in ; 

a 

where a is the planet's semi-major orbital axis. 

Although there is some evidence for a "pile-up" of planets at ~3 day periods (Cumming 
et al. 2008), the intrinsic distribution of planets with periods between 1-10 days is not 
well understood. Therefore it is necessary to make an assumption in assigning the period 
distribution in the Monte Carlo simulation. Hartman et al. (2009) use a uniform logarithmic 
distribution. We test this distribution, as well as a uniform linear distribution in period. 

With the parameters M*, M P) R+, R P , cosi, and P assigned to each star, a transit 
depth and duration was calculated and a box-shaped transit inserted into each synthetic 
lightcurve. For each parameter set, we simulated 10 randomly-selected phases. We adopted 
a simple box-shape model as limb darkening effects have been shown to be negligible in this 
signal-to- noise regime (Gould et al. 2006). The transits were inserted into lightcurves that 
had already been detrended using Sys-Rem, whereas a transit signal in our survey would 
need to pass through the SyS-Rem detrending. We therefore carefully tested and optimised 
the number of Sys-Rem iterations we used so that we would not degrade a real astrophysical 
signal. Only those trends that are present in many light curves with the same temporal 
timestamps are removed. 



(4) 
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6.4. Detection Efficiency and Effective Number of Stars Probed 

The BLS algorithm (Kovacs et al. 2002) was applied to each synthetic lightcurve, using 
the same parameters used for the actual SuperLupus dataset (see Section 4). In total, 800,000 
lightcurves were generated for the Monte Carlo and searched for transit signatures. Transits 
were deemed to be detected if the selection criteria used in the actual SuperLupus search 
(see Section 4) were satisfied. 

The detection efficiency, e, is defined as the fraction of lightcurves in which the inserted 
transit signal is detected in the Monte Carlo simulation: 



where 5i is a delta function representing whether the transit of star i is detected (Si = 1) 
or not detected (^ = 0), and N sim is the number of simulated stars in the Monte Carlo 
simulation. As discussed in Section 6.1, 24% of our survey stars are giants (log(g) < 4.0), 
and were automatically assigned a Si = regardless of the other system parameters. We 
also note that only transiting systems were simulated, so e does not include the geometric 
probability of transit. 

To determine the frequency of hot Jupiters in the field of the SuperLupus Survey we need 
to determine how many dwarf stars we have effectively probed for transiting planets, N pr . 
This number depends on the geometric probability of transit, the number of stars monitored 
to a given precision (for SuperLupus this is 20,465), as well as the detection efficiency. We 
calculated the geometric probability of transit, /f eo , for each system generated in our Monte 
Carlo simulation. N pr is therefore given as: 



Both the detection efficiency and the effective number of dwarf stars probed are set out 
in Table 4 for a variety of period ranges and assumed period distributions of hot Jupiters. 
Also tabulated are the mean transit depth, transit duration, and mean transit probability for 
each set of Monte Carlo assumptions about the size and period distribution of hot Jupiters. 
The results are presented for all periods simulated (1-10 days), as well as for period ranges 
of 1-3 days, 3-5 days, and 5-10 days in order to allow for direct comparison with the results 
of Gould et al. (2006) and Hartman et al. (2009). 

In Figure 5, N pr is plotted as a function of planet radii for each of the four planetary 
radii simulated. As expected, more stars are probed in all period ranges for large radii 




(5) 




(6) 
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planets, and more stars are probed for shorter period planets. The difference between the 
uniform and uniform logarithmic period distributions is most marked for the 1-3 d period 
planets. 




0.8 1 1.2 1.4 

Planet Radius (Rj) 

Fig. 5. — The effective number of dwarf stars probed (N pr ) for transiting hot Jupiters in the 
SuperLupus Survey as a function of planet radius. The solid lines represent the assumption 
of a uniform period distribution, while the dashed lines represent the assumption of a log 
uniform distribution. The square symbols are for 1-3 d periods, the circles for 3-5 d periods, 
and the triangles for 5-10 d periods. 
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7. Frequency of Planets in the Survey 

The frequency of dwarf stars that host hot Jupiter planets, /, can be calculated simply 
using the Monte Carlo results and the SuperLupus Survey results: 

pv 

where n is the number of planets detected in the survey. 

With only one (or possibly two) planets detected, Poisson statistics are used to determine 
confidence intervals, or upper limits in the cases of no planet detections. The probability of 
detecting n planets when A planets are expected is: 

P(n,X) = e x -. (8) 

Since the underlying distribution of planets is unknown, a Bayesian approach is used, with 
a flat prior for the uniform distribution, and a log prior for the logarithmic distribution. 
The upper and lower 90% confidence intervals {(Thigh and <Ji ow , respectively) are calculated 
by solving: 

r hiah e -^d\ 

f~ ^dX ' 1 j 

JO n! 

The values of Ohigh and o\ om were determined numerically from Eq. 9, with a symmetric 
confidence interval such that 5% of the probability distribution was below a\ ow and 5% 
above Ohigh- The limits Ohigh and ai ow were then used to calculate the upper and low limits 
for 90% confidence intervals. 

Where no planets were detected in a period range, a one-sided 95% confidence upper 
limit is used so that it can be directly compared with the two-sided 90% limits. In that case, 
n — 0, and Eq. 9 simplifies to: 

1 _ e -ohigh = o.95 . (10) 

This method of calculating confidence limits is similar to that used by Gould et al. (2006) 
and Hartman et al. (2009). 

The follow-up work presented in Weldrake et al. (2008) reveals that one of the SuperLu- 
pus candidates, SL-07, is indeed a hot Jupiter planet. Seven other candidates can be ruled 
out as false positives, while SL-04 and SL-06 remain as candidates, although follow-up would 
be difficult due to their faintness. We therefore present our frequency results for the case 
where Lupus-TR-3b is the only transiting planet detected (n — 1) and also for the case where 
either SL-04 or SL-06 is also a transiting planet (n = 2). The results are set out in Table 5, 
again for both uniform and uniform logrithmic period distributions, and are tabulated in 
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the same period ranges as used in Table 4. Frequencies for each of the four simulated planet 
radii are calculated, along with the mean over these radii, which is Rp = 1.1. 

7.1. Comparison with Other Surveys 

There are four other transit studies of non-cluster fields that report statistics for hot 
Jupiter frequency. 

The Kepler Mission (Borucki et al. 1997) yields very high precision, near continiuos 
photometry for around 150,000 target stars in a selected galactic field. Howard et al. (2011) 
calculate the the frequency of planets within 0.25AU of solar-type hosts stars selected from 
the Kepler target stars. For the period range of 1.2<P<10 d, and radii 0.7< Rj <1.4, very 
similiar to those considered in this work, the frequency of hot Jupiters is given as 0.37% 
(Figure 4 in Howard et al. 2011). This frequency is based only on "solar type" Kepler 
target stars, defined as those with temperatures of T e ff = 4100 — 6100 K and gravities of 
4.0 < logg < 4.9. While this gravity cut is essentially the same as is used in this work, the 
temperature cut means that a more limited sample of stars are being probed. Coupled with 
the fact that these stars are already drawn from a set of targets selected using "detectability 
metrics" that suggested they were most likely to give a detectable transit signal for terrestrial- 
size planets (Batalha et al. 2010), the frequency of 0.37% is higher than one might expect 
for ground based survey of random Galactic dwarfs in the field. Indeed the Kepler result is 
at the high end of the frequency of hot Jupiters determined from our SuperLupus survey, 
just consistent within our uncertainty. 

The MMT transit survey of the open cluster M37 (Hartman et al. 2009) included a 
significant number of non-cluster field stars, and these were analyzed separately to determine 
the hot Jupiter frequency for Galactic dwarf stars. Assuming a Jupiter-sized planet, and 
a logarithmic period distribution, that study found that the fraction of field stars hosting 
planets with periods of 1-3 d is <0.8%, and those with planets having periods of 3-5 d is 
<2.7%. These statistics are upper limits only, as there were no confirmed transiting planets. 

Analysis of the OGLE-III survey for transiting planets puts the frequency of hot Jupiters 
at 0.14% for 1-3 d period planets, and 0.31% for periods of 3-5 d (Gould et al. 2006). These 
statistics draw on survey fields in both the Galactic plane and the Galactic bulge, and were 
derived based on five transiting planets that had been confirmed from the OGLE-III survey. 
They rely on an estimation of the efficiency of the "by-eye" selection that was used by the 
OGLE team to select transiting planet candidates. 

The Hubble Space Telescope (HST) was used to undertake a survey for transiting planets 
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in the crowded bulge region of the Galaxy (Sahu et al. 2006). Sixteen candidates were dis- 
covered from a sample of 180,000 stars monitored, however only one of these (SWEEPS-11) 
was bright enough to be confirmed in the usual manner using radial velocity measurements. 
Assuming that all the candidates are indeed planets, it was reported that 0.42% of bulge 
stars more massive than 0.44 M Q have giant planets with periods up to 4.2 days (Sahu et al. 
2006). 

Our result for 1-3 d period planets, / < 0.13%, is consistent with the limits from the 
deep MMT survey of Hartman et al. (2009). It is lower than the OGLE results, but consistent 
within the lower half of their 90% confidence interval. For 3-5 d period planets, our result of 
/ = 0.15lo;ii is consistent with the upper limit from the deep MMT survey. It is lower than 
the OGLE results, but still within the 90% confidence intervals. For 5-10 d periods we are 
able to place limits on the frequency of hot Jupiters of / < 0.93%, assuming a logarithmic 
distribution in period and also that SL-04 and SL-06 do not host genuine hot Jupiters. If 
one of these candidates is a transiting hot Jupiter this frequency in this period range would 
be 0.44^34%, closer to the Kepler figure of 0.37% but still lower than the RV survey results. 
The comparisons are summarised in Table 6. 

The frequency of hot Jupiters can also be determined from radial velocity surveys. 
Cumming et al. (2008) report a frequency of 0.65 ±0.4% from the Keck Planet Search, while 
Naef et al. (2005) report a frequency of 0.7 ± 0.5% from the ELODIE planet search. These 
comparisons are also set out in Table 6. Both of these radial velocity results are above our 
upper 90% confidence interval (0.56%), although the confidence intervals of the two survey 
results overlap. 

There are, however, major differences between transit surveys and these radial velocity 
surveys that could lead to discrepancies in the frequency derived for hot Jupiters. Radial ve- 
locity surveys are sensitive to planetary mass, while transit surveys are sensitive to planetary 
radius. Also transit and radial velocity surveys usually monitor different stellar populations. 
Radial velocity surveys typically monitor carefully selected bright, sun-like stars. More re- 
cently, bright M dwarfs have also been monitored (Forveille et al. 2009). Transit searches 
monitor every star in the survey field that has the requisite signal-to-noise to yield high 
precision photometry. In a deep survey such as the SuperLupus Survey, this will translate 
to a stellar population with a sub-solar mean mass (O.9M in the case of the SuperLupus 
Survey). If hot Jupiters of a given period are less frequent around lower mass stars, then 
the frequency of hot Jupiters from deep transit surveys will be lower than radial velocity 
surveys. 
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8. Discussion 

The frequency of hot Jupiters in the Galaxy is an important quantity that will ultimately 
provide a constraint on models of planet formation and migration. 

It has been suggested that the frequency of hot Jupiters in globular cluster environments 
is lower than that of field stars, and that this may be due to crowding or low metallicity 
affecting planet formation, migration, or survival (Gilliland et al. 2000; Weldrake et al. 2005). 
Our results indicate, however, that there is little statistical disagreement between hot Jupiter 
frequencies in cluster and non-cluster environments, even if one of the remaining candidates 
turns out to be a genuine planet. This is consistent with the work presented in (Saders & 
Gaudi 2011) which concludes that there is no evidence to support that open clusters have 
a lower frequency of hot Jupiters. 

Initial estimates for planet yields from transit surveys turned out to be far in excess 
of the actual discovery rate (Home 2003). One of the many factors that led to this over- 
estimation was the adoption of the hot Jupiter frequency derived from early radial velocity 
surveys, which as we have shown is higher than is found from transit surveys. The Kepler 
result of 0.37% (Howard et al. 2011) obviously provides a robust statistic for hot Jupiter 
frequencies due to the high detection efficiency of that survey. However the result should 
be adopted cautiously when calculating expected yields from typical transit surveys, as it 
is based on a sample of "solar type" stars drawn from Kepler target stars, rather than the 
ensemble field stars monitored by most blind transit surveys. 

The authors thank David Weldrake and Brandon Tingley, who initiated the original 
Lupus Survey with PDS. We also thank Grant Kennedy and Tom Evans for assisting in 
gathering data for the SuperLupus project. 

Facilities: SSO:lm (WFI) 
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Table 1. Candidate selection criteria 



Criteria Threshold 



BLS Signal Detection Efficiency SDE > 7.0 

Effective S/N of transit a > 10 

BLS-Period P £ 1.02 d 

P ^ 2.00 ±0.02 d 
Magnitude V < 18.8 

Single-event fraction frac. < 0.9 

Not a variable star Weldrake & Bayliss (2008) 

catalog 



Table 2. Properties of Initial SuperLupus Candidates 
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Table 3. Candidate Test Results 



ID 


Discriminant 


Conclusion 


SL-01 


Saturated close neighbor 


Statistical False Positive 


SL-02 


High r] p 


EcB 


SL-03 


Additional Photometric Analysis 


Blended EcB 


SL-04 


Transit depth (0.042mag) suggests EcB 


Candidate 


SL-05 


Additional Photometric Analysis 


EcB 


SL-06 


Low a (10.2) 


Candidate 


SL-07 


Lupus-TR-3 


Planet 


SL-08 


Transit correlated to Sj 


Statistical False Positive 


SL-09 


Transit correlated to Sj 


Statistical False Positive 


SL-10 


Transit correlated to Sj 


Statistical False Positive 
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Table 4. Monte Carlo results 
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Table 4 — Continued 



Period 


Planet 


Mean 


Mean 


Mean 


Detection 


Stars 


Range, 


Radius 


Depth 


Duration 


probability 


Efficiency 


Probed 


Distribution 


(Rj) 


(mag) 


(hr) 


of transit 


(e) 


(N pr ) 




1.0 


0.016 


2.79 


0.066 


0.14 
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1.2 


0.023 


2.83 


0.068 


0.26 


281 




1.4 


0.031 


2.92 


0.069 


0.38 


431 



a Uni=Uniform distribution, Log=Logarithmic distribution 



Table 5. Frequency of hot Jup iters 



Period Planets Planet Frequency (/) 

Range & Detected 15 



Distribution a 


(n) 
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[0.20l°; 7 3 7 %] 
<0.23% 
<0.57% 
<2.13% 
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<0.39% 
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<1.09% 


<0.31% 
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<0.77% 
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(0.12^%) 
<0.15% 
<0.43% 
<1.41% 

[0.6615;™%] 


<0.18% 

<0.11% 
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<0.76% 


<0.14% 
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<0.50% 


o.iol^ 8 7 % 

(0.20l°f%) 
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a Uni=Uniform distribution, Log= Logarithmic distribution 

b Statistics in square brackets are for case where either SL-04 or SL-06 is a planet 



Table 6. Comparison of Planet Frequencies 



Survey & 
Reference 


Method & 
Population 


Assumptions 


Planet 
Frequency 


Kepler 

(Howard et al. 2011) 


Transit, 
"solar-type" 
dwarf stars 


1.2 < P < 10 day 
0.7 < R P < lARj 
n = 22 


0.37% 


Keck Planet Search 
(Cumming et al. 2008) 


Radial velocity, 
nearby dwarf stars 


M p > 0.5 Mj, 
P < 5 day 


0.65 ± 0.4% 


ELODIE Planet Search 
(Naef et al. 2005) 


Radial velocity, 
nearby dwarf stars 


M p > 0.5 Mj, 
P < 5 day 


0.7 ±0.5% 


SWEEPS 

(Sahu et al. 2006) 


Transit, 
bulge stars, 
M* > 0.5 M 


P < 4.2 day 
n = 16 


0.41^ % 


OGLE-III 

(Gould et al. 2006) 


Transit, 

bulge and Galactic 
disk dwarfs stars 


n=5, 1.0 <R P < 1.25 Rj 
log P distribution, 
1 < P < 3 day 
3 < P < 5 day 


o.i4l8:Si % 

0.31^% 



Table 6 — Continued 



Survey & 


Method & 


Assumptions 


Planet 


Reference 


Population 




Frequency 



Deep MMT Transit, n=0, R = 1.0 Rj 

(Hartman et al. 2009) Galactic disk log P distribution 

dwarf stars 0.4 < P < 1 day < 0.3% 

1 < P < 3 day < 0.8% 

3 < P < 5 day < 2.7% 



SuperLupus 
(this work) 



Transit, 
Galactic disk 
dwarf stars 



n=l, R= 1.1 Rj 
log P distribution 
n=l (Lupus-TR-3b) 

1 < P < 3 day 
3 < P < 5 day 
5 < P < 10 day 
1 < P < 10 day 



<0.15% 
0.151°;^% 
<0.93 % 



0.10 



-0.27 
-0.08 70 



% 
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